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00 We discuss the general link between mode-coupling like equations (which serve as the 
basis of some recent theories of supercooled liquids) and the dynamical equations governing 

CN mean-held spin-glass models, or the dynamics of a particle in a random potential. The 

^ physical consequences of this interrelation are underlined. It suggests to extend the mode- 

t-h coupling approximation to temperatures well below the freezing temperature, in which aging 

^ effects become important. In this regime we suggest some new experiments in order to test a 
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the short and long (a) time regimes. 



LPTENS preprint 95/47 S ACL AY preprint S95/068 



Typeset using REVTgX 
*Unite propre du CNRS, associee a l'Ecole Normale Superieure et a l'Universite de Paris Sud 



I. INTRODUCTION 



Let us face it: there are not so many techniques to deal with the score of strongly non- 
linear problems that Nature perversly offers, to the theoretical physicist dismay. Among 
others, one may of course cite fully developed turbulence 1 , but also interface growth and 
disordered systems 2 and strongly interacting liquids (i.e. glasses) 3 . The core of most of these 
problems is a non-linear dynamical equation, which we write in a symbolic way as: 

= -p(*)#M) - 9F(cf>) + 77 (I.f ) 

where <j)(x, t) is a vector field, F(<j>) is a non-linear (though generally local it can also be 
non-local) coupling term and rj a Gaussian white noise. The term containing fi(t) is a 
restoring force. We leave open the possibility that it can become time dependent in order 
to include in our study the cases where one imposes a 'spherical' constraint on the field </>, 
such as 4>(t) ■ 4>(t) = f. The coupling constant g serves as a book-keeping parameter to 
set up a perturbative expansion. This expansion can either be well-behaved or ill-behaved 
depending - say - on the dimension of space. It is in any case rather useless when g is of 
order I if it cannot be resummed in one way or another. A very popular and versatile class 
of resummation schemes amounts to performing a 'one-loop' self-consistent perturbation 
theory. Depending on the context, self-consistent approximations of this type have received 
the names of 'Mode Coupling Approximation' (MCA) 4 for critical dynamics or liquids, or 
'Direct Interaction Approximation' (DIA) 1 for turbulent flows; to some extent the Hartree 
approximation also falls in this category, as well as the refined version called 'Self Consistent 
Screening Approximation' (SCSA) 7 . In the MCA-DIA for the problem described by the 
Langevin process (I.I) one expands the relevant physical quantities to lowest non-trivial 
order in g and then replaces the bare objects in the correction term by the fully 'renormalised' 
objects that one wishes to compute. This amounts to resumming a particular (infinite) set 
of terms in the perturbation expansion. In this way, non-trivial self-consistent equations are 
obtained, which enable one to peep into the strong coupling regime. 

The problem is of course to try to control this procedure. An important step in this 
direction is to identify a model for which the self-consistent equations are exact (just as the 
Hartree approximation describes exactly the large N limit of an N component field problem). 
This is interesting for three reasons: first of all, it shows that if the underlying model is well 
behaved, the approximation does not violate any physical constraint. Second, the ingredients 
needed to build the model shed light into the physical content of the approximation. Third, 
one may hope to find a systematic expansion around this approximation. One can discuss 
in particular whether the interesting features of the self-consistent equations are or not an 
artifact of the approximation itself. 

This general concern is particularly relevant within the context of supercooled liquids, for 
which the Mode-Coupling Theory (MCT) t offers (at present) the most comprehensive and 



1"The Mode- Coupling Theory of glasses takes as a starting point an exact Liouvillian description of 
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successful description 3 ' 5 . It was understood long ago by Kraichnan 6 that the DIA approxima- 
tion for turbulence becomes exact when one considers a generalisation of the Navier- Stokes 
equation which contains some quenched disorder. Recently, it has been understood that this 
same approximation also becomes exact for a system with deterministic, but highly chaotic 
interactions 8 , which in fact are not very different from random ones (we shall return to 
this paper later on). This also holds for the simplest mode coupling equations with cubic 
interactions 9 . The existence of an underlying disordered problem is in fact a very general 
result: we shall show below that the MCA for a general non-linear F(<j>) and the dynamical 
generalisation of the SCSA 7 are the exact equations describing a suitably chosen disordered 
system. 

One extra difficulty of modelling 'true' glasses (with respect to spin-glasses) is that the 
effective disordered potential slowing down the particles is 'self-induced' by the dynamics 
itself, rather than arising from an external source of quenched randomness. At the same time, 
glasses and spin-glasses behave very much in the same way, suggesting that the difference 
between 'self-induced' and quenched disorder might not be crucial, at least in a restricted 
time window. This scenario has been substantiated within several mean-field like models 
in the recent years 10 12 . In a sense, the MCT introduces some quenched randomness into 
the glass problem, without specifying it explicitly. This might be a clue to understand the 
success of the MCT. 

The fact that MCT equations become exact for some disordered system suggests how 
to extend it to low temperatures, i.e. inside the glass phase. The MCT for glasses usu- 
ally addresses the temperature regime above the glass transition, in the supercooled liquid 
phase, where the property of time translation invariance (TTI) holds. This means that the 
correlations between time t and t' depend only on the time difference t — t' (as a matter of 
fact, the MCT is generally formulated directly in frequency space). 

However, as is now well documented experimentally in the case of spin-glasses 13 and 
other structural glasses 14 , this property does not hold in general in the glass phase. There 
is a non- vanishing 'waiting time' dependence in the correlation and response function - the 
'aging' effect. It was recently observed 16 that even such simple disordered models (as the 
ones for which MCA is exact) have a low-temperature out of equilibrium dynamics that 
is both soluble and indeed captures aging phenomena in qualitative agreement with the 
experiments. Hence it is important to know in general how Mode- Coupling-like equations 
can be written in a two-time formalism, without assuming TTI, since this allows one to 
make predictions deep into the glass phase, and not just above it. 

Finally, one could hope that some sort of perturbative expansion, taking the disordered 
system as a starting point, would bring one back closer to the original model, in particular 
accounting for finite dimensionality effects. 

The aim of this paper is threefold. We first show that the MCA for a general F(<j>) is 



the interacting particles but not Langevin noise r]. Through a series of approximations, similar in 
spirit to, but different from, the MCA, one obtains the so-called Mode-Coupling equations discussed 
in section 4, which happen to be identical to the MCA equations deriving from Eq. (1.1). There is 
thus a slight distinction between MCA and MCT. 
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equivalent to studying a general spin-glass system (Section II). Second, we show (Section 
III) that Bray's SCSA for the usual </> 4 theory amounts to studying a disordered version 
of the Bernasconi model 21 , which was studied recently precisely to give some flesh to the 
idea of 'self-induced' disorder in glasses. Finally, we summarize in Section IV the known 
results 16-18 on these disordered models and rephrase them in the context of supercooled 
liquids. We suggest that well controlled aging experiments deep below the dynamical glass 
transition temperature might serve as a crucial test for the Mode-Coupling description of 
glasses. The reader who is interested in the physical aspects of the discussion and less by 
the technical details can jump directly to Section IV. 



II. MODE COUPLING APPROXIMATION AND DISORDERED SYSTEMS 



We hrst describe the MCA on the simple case of a single scalar degree of freedom </>, with 
an energy 

H = ^f^ + ^\ (II.1) 

We assume that the dynamics of cf> in contact with a heat bath is described by the Langevin 
equation: 

^ = -Kt)<f>-§if + v (n.2) 

with initial condition cf>it = 0) = 0. The thermal noise rj is a Gaussian noise rj with 
(r)(t)) = and ( rj(t) r)(t') ) = 2T 6{t — t') (in the following the brackets will always denote 
an average over the realisations of the Gaussian white noise rj). 




Figure 1. Diagrammatic representation of the perturbative solution to Eq.(II.l). 
Crosses indicate noise and oriented lines indicate the bare propagator Go 
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Setting Go = [(J<(t) + J7] 1 , the perturbative expansion for <f>(t) is easily written as: 

9 



<t>{t) = Go ® f] - |j Go ® {Go ® f] • Go ® ?y • Go ® f]} + ... (II.3) 

where ® means a time convolution: (Go <E> /)(£) = Jo dt'Go(t,t')f(t') and • is a simple 
product. For the specihc form of Go in Eq. (II. 2), one has Go(t,t') = exp f t , dr fJ,(r)j . 
Eq. (3) can be graphically represented as in Fig 1. Crosses indicate noise and oriented lines 
indicate the bare propagator Go- 

Two quantites of interest are the (two-times) correlation function C(t, t') and the response 
function G(t,t') defined as 

C(t,t') = (<f>(t)<f>(t')), (II.4) 

G(M0^(f^> = ^<<^h(O>- (n.5) 

where the last equality holds for a gaussian noise. The diagrammatic expansion of C, G is 
represented in Fig. 2. 



C(M') = X + 3 + 6 > >< < + 



JL 
3! 




Figure 2 - a. Diagrammatic representation of the perturbative expansion 
of the auto-correlation function C(t,t') = (<f)(t)<f)(t')) 



G(M') = > + 3^^ >— + 18 > X > + 



JL 
3! 




Figure 2 - b. Diagrammatic representation of the perturbative expansion 
of the response function G(t,t') = 1/(2T) ((f)(t)i](t')) 
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In what follows we shall assume that the mass is renormalised in such a way that all 
tadpoles (i.e. the second diagrams in Figs 2-a and 2-b) are already resummed. 

It is useful to introduce the kernels E(t,t') and D(t,t') through the Dyson equations: 

G(t, t') = G (t, t') + f dt x [ h dt 2 G (t, h) X(t u t 2 ) G(t 2 ,t') , (II.6) 
Jt' Jf 

C(M') = [' dh f dt 2 G{t,h) D(t u t 2 ) G(t',t 2 ) . (11.7) 
Jo Jo 





E oc 




Figure 3. Diagrammatic representation of the MCA. The hrst two lines represent the 
Dyson equations, Eqs(II.6)-(II.7), which define the kernels E and D. The last identity 
gives the value of these kernels within the MCA in the case of the <f> 4 theory. An 
oriented double line denotes the full response G, an oriented single line denotes 
the bare response Go, and a crossed double line denotes the full correlation C. 
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The MCA for this problem amounts to an approximation of the kernels E(t,t') and 
D(t, t') where one takes their values at order g 2 and substitutes in them the bare propagator 
Go and the bare correlation by their renormalised values. This gives the following self- 
consistent equations: 



which are represented in Fig. 3. This approximation neglects 'vertex renormalisation': It 
keeps for instance the diagram depicted in Fig. 4- a that represents a line correction, while 
leaving aside the diagram drawn in Fig. 4-b that represents a vertex correction. 

It will also be useful in the following to note that the Dyson equations can be recast, 
after multiplying by Gq 1 , into the form: 





G - 1 (g)G = J+ E®G, 
G- 1 (g)G = J D(g)G + E(g)C', 




where X is the identity operator. More explicitly, 



dG(t,t') 

dt 
dC(t,t') 

dt 



f i(t)G(t,t') + 6(t-t')+ I dt"X(t,t") G{t",t') 



fi(t) C{t,t') + / dt" D{t,t") G{t',t") + / dt" E(t,t") C{t",t') 







The delta-function imposes G(t,t ) = 1. 




Figure 4 - a. Exemple of a graph that is kept in the MCA 
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Figure 4 - b. Exemple of a graph that is neglected in the MCA 



The basic remark is that the diagrams retained by the MCA are precisely those which 
survive if one considers the following disordered problem. First, one upgrades cf> to an A— 
'colour' object where a = {1,2,..., A}. The equation of motion Eq. (II. 2) is generalized 
to: 

— j- = -jl(t) <f) a - \g ^ Jaf3 1 S <f>b<f>~/<f>6 + Va (11.13) 

/3<7<5 

with independent noises rj a . This equation derives from the Hamiltonian 

Hj = g J»p 1 s ^a^p^i^s (H-14) 

a</3<7<5 

The couplings 3 a {5~f& are independent Gaussian random variables of zero mean and variance 
Ja/3-yS = 1/A 3 ). In the large A limit, the correlation: 

1 N 

c M = jj E (MW)) (n.15) 

a = l 

(where the overline denotes the average over the random couplings Ja^s) and the response: 



G(t,,| 4|0 (ILl6) 

precisely obey the MCA equations, Eqs. (11.7,11.8,11.6). (It is important to notice that 
the random couplings are quenched, i.e. time-independent random variables.) The fact 
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that MCA equations are recovered can be seen either directly on the perturbation theory, or 
through the use of functional methods given in Appendix A. A simple physical interpretation 
can be obtained through the cavity method 19 ' 17 where one shows that, in the large N limit, 
any one of the <f> a evolves through an effective linear Langevin equation: 

-£L = -^(t) S a + / £(*, f) <j> a (t') + Ut) + V a (t) (11.17) 
at Jo 

where £ a (t) is an effective (Gaussian) noise, with correlations self-consistently given by 
(Ca(t)Ca(t')) = D(t,t'). This result, derived in detail in Appendix A, gives back the MCA 
equations for the two point functions. It also provides a precise recipe to calculate higher 
order correlation functions within the MCA. 

The first to notice that the MCA (for the case of a 'quadratic' dynamical equation) 
corresponds to the exact dynamical equations of a disordered problem with a large number 
of components was Kraichnan 6 in the context of the Navier-Stokes equation. The important 
property of the random couplings which is used in the derivation is that the couplings are 
independent Gaussian variables. 

In the case of J's with three indices, this can also be implemented using a deterministic 
construction of the J a f5~ n in terms of Clebsch-Gordan coefficients of an 0(3) symmetry 
group +. This was first noticed by Amit and Roginsky 22 , and has been recently extended for 
dynamical problems, 8 ' 23 ' 9 . 

Interestingly enough, this equivalence between MCA and a disordered system extends to 
an arbitrary non-linear coupling F(cf>) (see Eq. (I)). Expanding F(cf>) in power series 

OO El 

W) = £-i> r ( IL18 ) 

r = 2 V - 

the natural generalisation of the MCA (i.e. neglecting all vertex renormalisation) reads: 



oo jp2 
I r 

r = 2 



^ = 9 2 E 77^? ^ t ')] r " 1 G(t, f) , (11.19) 



OO El 2 

D(t, f) = 2T6(t -t')+g 2 J2^T [C(t, f 0] r • (H.20) 



r = 2 r - 



(Note that for r odd, there appears an additional 'tadpole' contribution in Eq. (11.19), 
which we have assumed again, that it has been reabsorbed into the mass term fi(t)). The 
dynamical equations within the MCA for this extended model are readily obtained inserting 
these expressions for E and D in Eqs. (11.11) and (11.12). 

These equations can again be obtained as the exact solution of a problem with quenched 
randomness, the problem of N continuous spins <f> a interacting through the Hamiltonian: 



+The behaviour of these coefficients as a function of their indices looks however extremely 'chaotic' 
- the difference between determinism and randomness is thus probably very thin. See the discussion 
in Ref. [11]. 
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H J [<f>]= 9j2 F r J2 J ai ...a r+1 <f> ai ---<f>a r+1 (11.21) 

r>2 a\ <...<a r+1 



and the Langevin equation 

d<f> a 8Hj[4>] 

"ST = " -*£T + (IL22) 

where J ai ,..a r+1 are quenched symmetrical and otherwise independent Gaussian variables 
normalized as: 

1 



K Ja 1 ,..a r+1 ) 2 — Ar „- (11.23) 



Therefore the mode coupling equations corresponding to an arbitrary nonlinearity F(cf>) 
describe exactly a spin-glass problem with arbitrary multispin interactions. Let us note that 
in order to be well defined, the model defined by the Hamiltonian Hj must be supplemented 
by a constraint preventing the field <f> a from exploding in an unstable direction set by the 
coupling tensor J ai ...a r+1 - A convenient constraint is 

^E0) = C(M) = 1 (H.24) 

JV a = l 

which can be implemented dynamically through a Lagrange multiplier, acting as a time- 
dependent mass fi(t) which must be self-consistently determined. Another possible regu- 

i 2 r-\- 1 

larisation is to add to Hj a term Ng' /2(J2 a ) 2 ■> with g' large enough. As a matter of 
fact, this term precisely generates, for r odd, the tadpole contribution in the expansion of 
the original cj) r+1 model, provided one chooses g' = g/ (2~ ( L ^-)\). Surprisingly, it can be 
checked that this value of g' is not large enough to suppress the instability of the disordered 
model. We are thus led to conclude that the plain MCA approximation (i.e. without im- 
posing an extra constraint) for - say - the </> 4 model leads to spurious instabilities, at least 
at low temperatures. A similar conclusion was reached in Ref. [26]. 

Interestingly enough, the disordered multispin Hamiltonian can also be seen 29 ' 17 ' 18 as 
describing a particle evolving in an N dimensional space in a quenched random potential 
Hj[(j)] = V[(f>] This random potential has a Gaussian distribution with zero mean and 
variance^: 

with 



oo 



= 2 



F 2 

r 

(r + 1) 



V(x) = g 2 J2 i.. , r iM x ' r+1 ■ ( IL26 ) 



^Note that the sign of V differs from the convention adopted in Ref. [18] 
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The general mode coupling equations (II. 11), (II. 12), (II. 19), (II. 20), are thus also the exact 
dynamical equations for the problem of a particle in a random potential in large dimension 
N. In this last context, they are often written 17 ' 18 in a differential form obtained after 
applying the operator Gq 1 : 

= -p(t) C(t, t') + 2T G(t\ t) + 

+ f dt" V'[C(t, t")] G(t', t") + f dt" G(t, t") V"[C(t, t")] C(t", t') (11.27) 



dG gf^ = + s (t ~ t') + f dt" G(t, t") V"[C(t, t")] G(t", t') (11.28) 

The physical consequences of this general equivalence will be fully discussed in Section 
IV. The extension of the mapping to a space dependent <f)(x,t) (or to a multicomponent 
field) is straightforward. Several interesting physical examples involve an equation of the 
type: 

86(k t) 00 p A 

n ' ' =-(uk 2 + ^(k,t)-j2 E ^A(fc|fci, k r )^k 1 ,t)...4(k r ,t) + i 1 (k,t) 

dt r = 2k u ..k r rl 

(11.29) 

where <p(k } t) is the Fourier transform of <f)(x,t), and rj(k } t) a Gaussian noise such that 
(r](k,t)r](k' ,t')) = 2T(k)S(k + k')S(t — t 1 ). The case of the KPZ equation 2 corresponds to 
r = 2, C 2 (k\ki } k 2 ) = [k\ • k 2 ] S(ki + k 2 + fc), while domain coarsening in the </> 4 theory 
corresponds to r = 3, £ 3 (fc|fci, k 2} k 3 ) = S(ki + k 2 + k 3 + with a negative /i 24 . The 
Navier-Stokes equation is similar to the KPZ case, with however an extra tensorial structure 
due to the vector character of the velocity field 1 . 

The correlation and response functions now become k dependent: 

S(k + k')C(k, t, f) = (]>(k, t)4>(k\ t')) (11.30) 
S(k + k')G(k,t,t') = (^^) (11.31) 

The generalized MCA equations then read (assuming that the structure factors 
) are invariant under the permutation of k r ): 

00 p2 

E(k,t,t') = g 2 Ez~T~[\! ^ C r (k\k u k r )C r (k r \k u k) 

r = 2 V /' ki,..k r 

C{k 1 ,t, *')...C(!b r _i, t, t')G{k r ,t, f) (11.32) 

00 p2 

D(k,t,t') = 2T(k)S(t-t')+g 2 J2^T E K)f 

r=2 ' ki,..k r 

C(k 1 ,t,t')...C(k r ,t,t') (11.33) 
where E(fc,t,t') and D(k,t,t') are defined in analogy with Eqs. (II. 6), (II. 7). 
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III. SCSA AND DISORDERED SYSTEMS 



Another useful resummation scheme is the 'Self- Consistent Screening Approximation' 
introduced by Bray in the context of the static <f> 4 theory 7 . It amounts to using an n 
component vector field <f> and resumming self consistently all the diagrams appearing in 
the large n expansion (n is the number of components of <f> n ), including those of order ^. 
This approximation can also be seen as a MCA when one introduces an auxiliary field by 
rewriting the Langevin equation for the </> 4 theory as: 

^ = -,■<*(( )-9<i(i )<r(i ) + ,,, (III.l) 

*<0 = *f (HI.2) 

with g = ly. (The factor 2 has been introduced for later convenience). In this form one gets 
back a problem similar to the ones studied before, which can be seen as two coupled fields 
<f> and a evolving with bare evolution operators 

(G^y 1 = //(*) + j f ■ (G a0 )=l (III.3) 

Once this fictitious decomposition of the non-linear coupling is performed, one can apply the 
MCA on the coupled equations (III. 2). Of course, if MCA were exact, the approximation 
would give the same results as in the previous paragraph. The fact that it is only approx- 
imate leaves room to a certain freedom on the starting point to improve (or deteriorate) 
the quality of the approximation (see Ref. [26] for a related discussion). Introducing two 
correlation functions C^{t,t') and C CT (t,t'), and two response functions G^, G>, together 
with the corresponding kernels E^, D^, E^ and D a , defined (separately for each field <f> or 
a) as: 

G' 1 = (Go)" 1 - £ ; C = G®D®G T (III.4) 

one finds the following result for the kernels**: 

E^(M') =g 2 C a {t,f)G 4> {t, t') - gC^(t } t')G«(t } t') (III.5) 
D^t, f) = 2T6{t - t') + g'C^t, t')C a {t, f) (III.6) 
E ff (M') = -gC+{t,t')G+{t,t') (III.7) 

D a (yt')= l -[C 4> (t,t')] 2 (III.8) 

It turns out that, again, these dynamical equations are exact for a certain (mean-field 
like) spin-glass model. Let us define the following 'spin-glass' Hamiltonian: 



**We have again reabsorbed the tadpole contribution in ji(t) 
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j N N 

« = ^£E^V»<M 2 + f££ (ni.9) 

ZJV A=l a</3 L » 

where the j£ ^ are identically distributed independent (apart from a constraint of symmetry 
in the two indices a, f3) random variables, such that j£ ^ equals 1 with probability and zero 
otherwise. This model was proposed and studied in Ref. [11] , as a disordered proxy of the 
Bernasconi model (which serves as a model for glassy behaviour without randomness). The 
calculations showing that the dynamics of the model defined by Eq. (III. 9), with J = g, 
exactly reproduces the self-consistent dynamical equations (III. 8) are given in Appendix 
B. The relation with Bray's SCSA can be directly seen on the statics of the disordered 
Hamiltonian, Eq. (III. 9). It is straightforward to show that, in equilibrium, 

CJt,t) = (111.10) 

K } 1 + #[CV(M)] 2 1 ^ 



: (III.ll) 

which indeed coincide with Bray's equations in zero dimensions with the identification n = 1, 
and his choice for ^ = 2. 

The tadpole term in the expansion of the </> 4 theory can also be taken care of by adding to 
the disordered system's Hamiltonian (III. 9) a term in g' /2N(J2 a ^a) 2 - We notice that within 
this approximation the energy (III. 9) is always positive, which ensures that the dynamical 
version of the SCSA are well defined, at variance with the MCA (cf. above). If the quadratic 
term in the original hamiltonian is positive, then the spin-glass sytem is unfrustrated: it 
has a single ground state at <f> a = 0. If instead we consider a double well </> 4 theory with 
a negative //, we find a frustrated spin-glass system. The usual Bernasconi model involves 
Ising spins with the same coupling as in (III. 9). It is recovered here in the limit where 
fi — > — oo. The study of the physical content of this dynamical SCSA is left for future 
work 25 . 



CJt,t) 



IV. PHYSICAL DISCUSSION: MODE COUPLING BELOW T G 



We have shown in the previous sections that the Mode- Coupling Approximation (or 
the SCSA) for a non-linear dynamical Langevin equation amounts to studying an auxiliary 
Langevin process for a system with quenched disorder. In particular, the MCA for the 
Langevin process (1.1) described by the non linearity F(cf>) leads to the pair of coupled 
dynamical equations for the correlation C(t,t') and the response G(t,t') written in (11.28). 
As we have seen, these equations describe exactly the dynamics of a particle in a random 
potential in a large dimensional space, or else as a certain type of mean-field spin-glass 
system with multispin couplings. 



12 



Actually the usual mode coupling equations which have been used successfully in the 
study of supercooled liquids are a special case of these general equations. The MCT is 
written in terms of the density-density correlation function which is normalized to one at 
equal times, i.e. C(t,t) = 1, corresponding to the spherical constraint (11.24). In the case of 
a supercooled liquid one studies a system in its high temperature phase where it obeys time 
translation invariance (TTI), together with the Fluctuation-Dissipation Theorem (FDT). 
The first of these properties allows to write the correlation and response as functions of time 
differences only: C(t,t') = C(t - f) and G(t,t') = G(t - f). The FDT, which states that 

G(r) = -| e(T)d T C(r) , (IV.l) 
(where r = t — t') } enables one to rewrite the mode coupling equations as a single equation: 

d T C(r) = -AooC(r) + i f dr"V[C(T-T")]d T „C(r") , (IV.2) 

1 Jo 

where fi^ = lim^^ ji(t) - 1/T V'(l). 

Eq. (IV.2) is basically the general Mode-Coupling equation for the density correlations in 
a supercooled liquid above the dynamical transition temperature introduced by Leutheusser, 
Gotze and others 3 as a 'schematic' model for the ideal glass transition. This similarity (in 
the high temperature phase) was already pointed out in 27 ' 28 . The only difference lies in 
the fact that the Mode-Coupling equations also possess an 'inertial' term <9 2 C(r). In the 
notation of Ref. [3], the 1 Ff models correspond to the case where the non linearity is a 
pure power law, where only F r is non zero; the F, n , r2 models correspond to a nonlinearity 
which is a sum of two powers with F ri , F r2 different from zero, etc. One should note that the 
Mode-Coupling equations for supercooled liquids were written from the start within a TTI 
formalism ^ . The analysis of these mode coupling equations (IV.2) for supercooled liquids 
has shown the existence of a dynamical phase transition at a certain temperature Td (which 
is traditionally called T c in the MCT), and identified two classes of behaviours (called A and 
B) when the temperature decreases and approaches Td- This same classification has also 
been discussed in the spin-glass dynamics framework, where the temperatures lower than 
Td (i.e. inside the spin-glass phase) has also been discussed and has led recently to several 
interesting developments. As noted recently in Ref. [9], these studies of spin-glass dynamics 
below Td provide a natural generalisation of the mode coupling equations below the glass 
temperature, the physical content of which we shall discuss. 

Let us thus summarize the important results associated to the dynamical equations 
(11.28) and rephrase them in the context of the Mode-Coupling theory. 

There exists a critical temperature Td (or a set of coupling constants F r ) separating a 
'liquid' (or paramagnetic) phase where TTI and FDT hold. The dynamics is described by 
Eq. (IV.2) and the correlations decay to zero at large times: C(r) — > when r — > oo. The 



t^If one attempts to extend directly (IV.2) to the low temperature phase keeping TTI and FDT 
as in Ref. [3], one obtains a theory yielding different predictions, the meaning of which is not clear 
to us. 
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transition can be of two types. In a first class of systems the transition is a continuous one: 
the analysis of the static situation through the replica method leads to a 'continuous replica 
symmetry breaking' 19 transition occuring at the temperature T s which coincides with the 
dynamical temperature Td where the ergodicity is broken. This corresponds to class A in 
the classification of Ref. [3]. The second class of systems have a very different behaviour 
where the static transition temperature T s is smaller than the dynamical one Td- This static 
transition, in the replica language, is a 'one step replica symmetry breaking' transition, 
which means that it is a first order transition from the point of view of the order parameter 
(but it is second order from the thermodynamic point of view). It corresponds to class B in 
the classification of Ref. [3]. We shall concentrate on this second category, which is supposed 
to be the most relevant for a description of the structural glass transition. In that respect, it 
is interesting to remark that class B systems correspond, in the equivalence with a particle 
in a random potential, to the case of short range correlations of random potential, whereas 
class A systems correspond to long range correlations 32 ' 17 ' 18 . 

Before describing the quantitative feature of the dynamical transition for class B systems, 
a few comments on their physical relevance is in order. The existence of a dynamic transi- 
tion above the static one is associated with the appearance of many metastable states and a 
breaking of ergodicity at T^, which does not reflect onto the equilibrium (Gibbs) measure 28 . 
However, this effect can exist only at the mean-field level, and it has been suggested that 
in finite dimensions some nucleation processes 10 ' 31 smooth the transition at Td and replace 
it by a crossover temperature range where the relaxation times will increase very fast with 
decreasing temperature. The glass transition temperature T g , empirically defined by the 
fact that the relaxation time (or the viscosity) reaches a certain conventional value, would 
therefore lie below the mean-field Td (but above the static transition temperature T s ). Ac- 
tually, the same type of argument has been developed in the study of supercooled liquids, 
where some 'activated processes' are supposed to smooth out the dynamical transition 3 . 

Hereafter we shall first recall the existing results for the dynamics above Td in spin-glasses 
and in supercooled liquids. These lead to the well-known predictions of the mode coupling 
theory for the relaxation just above Td, which have been tested experimentally. Then we 
shall recall the results of spin-glass dynamics below the dynamical transition. These lead 
to some predictions for the (off equilibrium) dynamics which should apply to glasses at 
much lower temperatures (smaller than T 3 ), such that the relaxation time is larger than the 
experimental time scale. 

- For T > Td } the analysis of Eq. (IV. 2) is sufficient. One finds that 29 , for T close to 
(but above) T c , C(r) has the form given in Fig. 5, with a plateau and the celebrated a and 
f3 regimes, characterized by two exponents a and f3 related through 3 : 

r 2 [i + «] r 2 [l-/3] T V"(q) 



(IV.3) 



r[l + 2a] r[l-2/3] 2 (V"(?)) 3 / 2 

with the value of the correlation at the plateau q given by (1 — q) 2 V"(q) = T 2 . Note that these 
two exponents a, f3 are usually called 6, a in the Mode-Coupling litterature. We however feel 
that it is more appropriate to call a the exponent corresponding to the a peak, and f3 the 
one corresponding to the f3 peak! 

- For T < Td } there appear diverging relaxation times in the problem. It has been realised 
recently 16 that in this case one needs to take into account carefully the existence of an initial 
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time for the dynamics. Stated differently, one must abandon TTI, as the age of the system 
becomes an important time scale in the problem. This leads to the existence of so called aging 
effects which have been observed in spin-glasses 13 , polymer glasses 14 and also in a variety 
of other systems 33 . A study of the full dynamical equations (11.28) shows that one must 
also abandon the FDT. The system is out of equilibrium, but one can nevertheless obtain 
some information on its behaviour. The correlation function C(t,t') (and similarly G(t,t')) 
must be decomposed into two parts (see Fig. 6): C(t w + T,t w ) = Cfdt(t) + C(t w + T } t w ) } 
Cfdt is TTI, it is related to Gfdt through the FDT, Eq.(IV.l), and corresponds to the 
high frequency dynamics (/3 peak), while the aging part C(t w + T,t w ) is a function of the 
ratio A = h(t w + r)/h(t w ) only. The 'effective time' h is still not determined theoretically, 
but a likely possibility, advocated in 15 , is that hit) = t. In other words, the relaxation 
time corresponding to the aging part of the correlation is the experimental waiting time 
t w itself. The a regime thus still exists for T < T c if the waiting time is finite. Only 
in the limit t w — > oo will the correlation relax to a non zero value. This is the 'weak 
ergodicity breaking' scenario proposed in Refs. [15], [16]: lim T _ ) . 00 linii m _ ) . 00 C(r + t w ,t w ) = q 
and lim T _ ) . 00 C(r + t w ,t w ) = 0, \/t w finite. The exponents a and f3 are thus still well defined 
for finite t w (see Fig 6) as 

C(t w + T,t w )^q + c P T- p if C^q (IV.4) 



T 



C(t w + T,t w )~q-c a ( — ) a if C^q. (IV.5) 

Aging is manifested in the t^-dependence of T w = dt w /d\n(h(t w )), which is an increasing 
function of t w . In the simple case where h(t w ) = t W} one has T w = t w . 

The exponents a and f3 are now given by a modified relation which reads 18 : 

i r 2 [i + q] _ r 2 [i-^] _r v"'(g) 

r[l + 2a] r[l - 2/3] 2 (V"( g ))3/2 v • ; 

with q given by (1 — q) 2 V"(q) = T 2 . x is a temperature dependent number, < x < 1. A 
crucial observation is the fact that this number is not arbitrary and could be in principle 
measured. It actually provides the quantitative measure for the violation of FDT. More 
precisely, x is defined as 16 : 

£(M0 = ^7^, (IV-7) 

where we assume t' < t. The usual FDT relation would state that x = 1. Glassy dynamics 
below T g gives a value x < 1, which also governs the relation between the exponents a and 
f3 in (IV.6). 

We shall not expand here on the case of class A situations, but just mention that the 
behaviour in the low temperature phase is more complicated 16-18 . The correlation and 
response have to decomposed into two parts as in class B situations but the behaviour of 
the aging parts C, Q cannot be characterised by a single function hit) and the violation of 
FDT is given, in the limit of large times, by a non-trivial function of the correlation function 
X[C] (instead of the single constant x). 
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T 

FIG. 5. Decay of the auto-correlation function above the critical temperature. 
C(t + t w ,t w ) = C(t) vs t. 

Let us finally say a word on the distinction between explicit and spontaneous non equilib- 
rium. Throughout this paper we have discussed extensions of mode-coupling-like dynamical 
equations which reduce to the usual ones if one assumes that TTI and FDT hold, as when the 
system they describe has achieved equilibration in some component of phase-space after some 
finite transient. However, we now know that such equations may admit a low-temperature 
glassy phase in which the equilibration time is infinite: there is violation of TTI and FDT 
at arbitrarily long times. The reason why this spontaneous non-equilibrium happens is that 
the equilibration time diverges, or at least becomes extremely large, with the system size. 
On the other hand, there are systems such as surface-growth (described by the KPZ equa- 
tion mentioned above) and stirred turbulence which are by construction non-equilibrium 
situations; their equations of motion do not admit any equilibrium solution even for a finite 
system. One can then wonder how to recognize if a given set of equations for response and 
correlation functions has explicit or spontaneous long-time non-equilibrium. It is interesting 
to notice that this question has a clear meaning within the supersymmetrical field theory 
described in Appendix A for the dynamics of disordered system. Any Langevin process 
that derives from a potential automatically yields an action that possesses a certain (su- 
per)symmetry (spontaneously broken if there is a glassy phase), while systems with explicit 
non-equilibrium have a dynamical action that break this symmetry explicitly. 
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T 

FIG. 6. Decay of the auto-correlation function below the critical temperature. C(r + t w ,t w ) 
vs. t for different waiting times, t wl < t w2 < t w3 . T w = dt w /dln(h(t w )). 

V. SUMMARY AND CONCLUSIONS 

Summarizing, the major prediction of the Mode- Coupling theory of glasses for their 
super-cooled liquid phase is the existence of a critical temperature Td below which the 
correlations do not decay to zero, and above which one observes two relaxation regimes (a 
and characterized by a power-law behaviour with exponents related by Eq. (IV. 3) - 
which is indeed qualitatively consistent with experimental data 3 ' 5 . However, a quantitative 
comparison is difficult since experimentally, the relaxation time r(T) does not diverge at 
Td but grows rapidly (a la Vogel-Fulcher) as the temperature is decreased further. In the 
Mode-Coupling approach, this is ascribed to some 'activated (or jump) processes' which 
must be taken into account in a phenomenological way. This can be rephrased differently: 
since we have argued that the Mode- Coupling equations are equivalent to the dynamics of 
a mean-field disordered model, it is to be expected that actual finite dimensional systems 
should depart from this ideal behaviour. A nucleation-like mechanism was proposed in Refs. 
[30], [31] to account for the smearing out of the transition in finite dimensions, but a detailed 
understanding of this mechanism is still lacking. This is in some sense related to the general 
question of assessing the quality of the MCA, and constructing perturbative schemes to 
move away from it 8 ' 20 . 

In order to by-pass this difficulty brought about by a finite relaxation time scale be- 
low Td, we propose that experiments should be done below T g , so that the experimental 
time scales t w are much smaller than t(T). Experimental protocoles should allow one to 
monitor in a systematic way aging effects (i.e., the fact that the correlation function does 
depend on t w itself), and to obtain the curves corresponding to Fig. 6. The crucial test of 
Mode-Coupling theory would then be to measure both correlations (or noise) and response 
functions (such as dielectric properties or elastic moduli) to observe the violation of FDT 



17 



and check Eqs. (IV. 6) and (IV. 7). It should be emphasized that most ol the experimental 
data on supercooled liquids (and spin-glasses for that matter) can alternatively be inter- 
preted within a phenomelogical model of 'traps' 34-36 ' 15 ' 37 ' 38 , where each particle diffuses in a 
random potential created by its neighbours. It would be interesting to understand the pre- 
cise relation between this phenomenological picture and Mode-Coupling equations 39 , which, 
as we have discussed, also describes a particle in a random potential, albeit in infinite di- 
mension. In any case, the genuine non-trivial prediction of the Mode-Coupling theory is 
that the equilibration process within a 'trap' (described by the exponent f3) and the aging 
process involving jump between traps (described, at least for small t/t w , by the exponent 
a) are intimately related through Eq. (IV. 6). This is why we believe that its investigation 
is worth the experimental effort. 
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Appendix A: 

Symmetries and the dynamical equations using the functional supersymmetric 

formalism 



In this Appendix we review the technique of superfield notation, which is useful for three 
purposes: 

• It provides a direct dictionary between statical and dynamical developments. There 
is, in this notation, a one-to-one correspondence between static and dynamical diagrams, so 
that one can for example talk indistinctly of 'static' and 'dynamical' MCT or SCSA. 

• It makes the diagrammatic developments simpler, grouping 2 L ordinary diagrams of L 
lines into a single superdiagram. 

• It makes explicit a supersymmetry (SUSY) of the action (or equations of motion) 
which embodies the equilibrium theorems. It is then possible to see directly from the form 
of the action whether there is explicit or (possibly) spontaneous non-equilibrium phenomena 
(breaking of SUSY), and to check that an approximation scheme or an effective theory does 
not spoil artificially the possibility of equilibrium. 

We start from a Langevin equation: 

d<f> a dH 

where rj a (t) are Gaussian random variables with zero mean and variance (r] a (t) rjp(t')) = 
2T 6 afj 6{t-t'). 

We now construct the Martin-Siggia-Rose 43 ' 44 functional for the expectation value of an 
operator 0{cf>) as: 



w»=/*]wn«(^^-uo) ^ 



d_ g d 2 H 

dt a,fj d<p a d<pp 



(A.2) 



where the measure D[<f>] is defined as D[<f>] = Yl a D[<f) a ]. 

Exponentiating the delta function through Lagrange multipliers <f) a (t) and the determi- 
nant through anticommuting variables ('ghosts') £ a (t) ,£ a (t), and averaging away the noise 
we obtain ++ 



++The precise meaning of Eq. (A. 4) is seen by going back to the Hilbert-space problem 42 of which 
it is a functional representation. Then, Eq. (A. 4) represents an imaginary-time evolution opera- 
tor associated with the Hamiltonian: J2 a Pa(Tp a — i§^r) + 2 a « d!j> 90 ^ a Pi wnere Pa = — ^gf - The 
original Fokker-Planck process is recovered by restricting the problem to the zero-ghost subspace, 
i.e. by considering diagrams without fermionic legs. The problem of which convention (Ito or 
Stratonovitch) is used is simply the usual problem of factor orderings in the functional represen- 
tation. Eq. (A. 4) with the assumption that 4> a (t)^-(t) is understood as (f) a (t + )j^-(t) is then an 
unambiguous representation of the Hilbert-space problem (in the Ito convention). 
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(0(<f>)) = / D[<f>] D[<f>] D[£] D[£] 0(<f>) exp(-S) 



(A.3) 



where 



s 



dt 



8H 8 2 H 



(A.4) 



The expression for S can be written in a compact form in superspace introducing two 
anticommuting Grassmann variables 6,6: 



[6,6]+ = 6 2 = 6 2 = . 
The integrals over these variables are defined as: 



U6 



U6 = 



6d6 



6d6 = 1 . 



(A.5) 



(A.6) 



The fields, Lagrange multipliers and ghosts are then encoded in the (bosonic) superfield: 



= Mi) + & Ut) + Ut) + <j> a {t) 6 .6 
Using Eqs. (A.5)-(A.6) and (A. 7) one obtains, in terms of the superfields $ c 



(o) = |n ^t $ «] ° ex p / 31 



(A.7) 



(A.8) 



where we have denoted 1 = (6 } 6 } t), dl = d6 d6 dt and the differential: 



8 2 8 2 
=2T-~ +26 



8 



86 86 



86 dt dt 



(A.9) 



The important point about expression (A.8) is that, apart for the first 'kinetic' term in the 
exponent and the integration over the 'time-like' coordinates dl = d6d6dt, the rest has the 
same form as the partition function. Furthermore, the correlation function between two 
superfields: ( $ a (f )$&(2) ) (with f = 2 = 6 2} 6 2} t 2 ) encodes all correlations and 

response functions. 

Consider now a single superfield $(f) = <f>(t) + 6£(t) + £(t)6 + cf>it)66. For a system 
satisfying causality, the non-zero expectation values of the autocorrelation function are: 

Q(f,2) =($(1)$(2)) 

= (Hti)Hh)) + (0 2 - h) [02 (<f>(t 1 )}(h)) + e 1 (Ht 2 )Uh))} 

= C(h, h) + (6 2 - 6 X ) [6 2 G(h,t 2 ) + 6 X G(t 2 ,h)] (A.f 0) 

Before going into diagrammatic computations, we need to define convolutions of two- 
point functions, as in: 



Q c (l,2) = Q a ® Q b 



d3 Q a (f,3)Q 6 (3,2) 



(A.ff) 
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and simple products, as in: 



g c (i,2) = g a (i,2).g 6 (i,2) 



(A.12) 



In what follows we shall denote with • any function based on usual products (e.g. Q' 3 = 
Q • Q • Q), and we shall omit ® when writing convolutions (e.g. Q 2 = Q (/) Q). 

At any step, one can go down to the 'components' of Q. For example, if Qi (i = a, b) are 
of the form (A. 10), 

Q l (l,2) = C l (t 1 ,t 2 ) + (9 2 -9 1 ) [0 2 G l (t 1 ,t 2 ) + 9 1 G l (t 2 ,t 1 )] (A.13) 
then Q c (l, 2) obtained with the two products above is of the same form with: 



(A.14) 



C c (t 1 ,t 2 ) = J dt'iCaih^Gbih^ + Gaih^Cbit',^)] 
G c (t 1 ,t 2 ) = I dt' G a {t u t')G h {t\t 2 ) 



for the convolution (A. 11), and: 



C c (ti,t 2 ) — C a (ti,t 2 )Cb(t 2 ,t 2 ) 

G c (ti,t 2 ) = C a (ti,t 2 )Gb(ti,t 2 ) + G a (ti,t 2 )Cb(ti,t 2 ) 



(A.15) 



for the usual product (A.12). 

The Mode Coupling approximation 

Let us now turn to the example of equation (II. 2). The Hamiltonian is given by: 

m = ^fV + (A.i6) 

The dynamical functional then reads: 

(0) = f D[$] 0($) exp -jdl + Kt)) + ^ 4 ) • (A-17) 

Diagrams are constructed as usual. They are based on the Gaussian integrals with bare 
propagator 



Go = [- J D( 2 ) + /i(t)]- 1 

defined by: 

Go ® + ^t)] = 6 

where we have defined the superspace delt 

S(l-2) = S(t 1 -t 2 )(9 1 - 9 2 )(9 1 - 9 2 ) . 



(A.18) 



(A.19) 



(A.20) 
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With these definitions, a 'superdiagram' is obtained just as an ordinary diagram, with 
now the labels 1,2 encoding both the times and the Grassmann coordinates, indicating 
the 'components' C,G of Q. Each line of the superdiagram stands for a line of field-field 
( </>(ti)</>(t 2 ) ) contraction (~ C) and a line of held-noise ( ^(^1)77(^2) ) contraction (~ G). The 
Dyson equations (II. 7)- (II.6) are both encoded in (see Fig. 7): 

Q = Go + G Cx)£®Q (A.21) 

The (super) mass-operator S can be now calculated within the mode-coupling approxima- 
tion, which again consists in neglecting vertex corrections. It is given by: 

E(1,2) = ^Q' 3 (1,2) (A.22) 

Introducing this into (A.21), and multiplying (in the sense of convolution) by G _1 , we 
obtain 

~d[ 2) Q(1, 2) = -fi(t) Q(l, 2) + 6(1 - 2) + ^ J d3 [Q(l, 3)]* 3 Q(3, 2) (A.23) 
This is the equation of motion. For a general non-linear coupling, one similarly obtains 

-D[ 2) Q(1, 2) = -fi(t)Q(l, 2) + 6(1-2) + f d3 V'*(Q(1, 3)) Q(3, 2) (A.24) 
which, in components, is nothing but Eqs. (11.27) and (11.28) . 



Q G G f ' ^ Q 

= = = = = = + H E 1= = = = 



// N\ 

// w 

\\ // 



\ I 
\ / 



Figure 7. Diagrammatic representation in SUSY formalism of the Dyson Equations. 
In this notation one encodes the C — G diagrams of Fig. 3 in only one super-diagram. 
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The disordered model for the MCA 



One can reconstruct an action functional of which (A. 24) is a stationary point. To do 
this, we multiply (A. 24) to the right by Q~ x : 



(D? - v(t))S(l - 2) + Q~\l, 2) + V"(Q)(1, 2) = 
Which can be written as: 

^ = 
6Q 

2S= f dl d2 \{-D { ? + //(*)) Q(l, 2) - V*(Q)1 - Tr Ln[Q] 



(A.25) 



(A.26) 



Let us now see how these are the exact equations of motion for a disordered system. The 
dynamical generating functional reads in superspace notation 



/ di [f;i$ a (i)(-^ 2) + M*))Mi) + ^[*] 



D[$] exp 
with the disordered Hamiltonian 

oo 

r>2 ai <-..<a r+1 

correlated as in (11.25). 

Averaging over the couplings (see Eq.(II.25)), we obtain 



D[$] exp 



exp 



Introducing the order parameter 



1 



Q(1,2) = ^EW.(2) 



1 



through 

1 = J D[Q] D[Q]exp 
yields 

' D[$]D[Q]D[Q] exp 

,(2) 



dl d2 NQ(1, 2)Q(1, 2) - Q(l, 2) £ $a(l)$a(2) 



J1J2 iVQ(l,2)Q(l,2) - V(Q(1,2)) 



E Ml) j + MO) *(1 " 2) + Q(l, 2)) $ a (2) 



(A.27) 



(A.28) 



(A.29) 
(A.30) 

(A.31) 



(A.32) 



(A.33) 
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We can now make the shift: 

Q = (-D? ) + f i(t))S(l-2) + Q(l } 2) 
and the integration over $, to obtain 



(A.34) 



D[Q]D[Q] exp 



TV 
exp 



dld2 (Q(l, 2)Q(1, 2) + (-D { ? + //(*)) Q(l, 2) - V*(Q(1, 2)) 

TV 



Tr Ln Q 



(A.35) 



(A.36) 



Using saddle-point evaluation, we can eliminate Q, and obtain: 
D[Q] exp(-NS) 
2S= J dld2 [(-D[ 2) + /,(t)) Q(l,2) -V'(Q(1,2))} - Tr Ln[Q] 

The saddle-point equation over Q yields (A. 25). 
SCSA and the Bernasconi model 



In a similar way we show that the equations of motion for the disordered Bernasconi 
model (III. 9) coincide with the equations arising from the SCSA applied to the </> 4 model. 
The generating function for the dynamics reads, in superspace notation: 



D[$] exp 



r ( 1 N j N / N \ 2 \ 

J y Z a=l ZJV A=l \a</3=l / y 



(A.37) 



Making a Gaussian transformation by means of superfields cry. 
£>[$] D[(r] exp 



„ / i N 

' a</3,A ^ ' 



(A.38) 



Averaging over the couplings: 

j D[$]D[(t] exp 

1 

exp 



1 



N 



a = l 



-\ Jdld2 (^2J (^E ^(1)^(2)) +*(l-2)j $>A(lM2)j 



(A.39) 



Note that putting 



< a/3 = 0(1 /N) creates an infinitely 

strong antiferromagnetic force. One can then neglect in the previous expression a term 
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J dld2 ^Ea ^a(l)/V^) ^/j(2)/V^V) since it is of order 0(1) and hence much smaller 

than the other 0(N) terms. One can just as well put p = and see directly how this 
term disappears. 

Introducing two order parameters Q§ and Q a as in (A.31)-(A.32) 



iVQ $ (l,2) =£$ a (l)$ a (2) NQ a (1,2) = £^(1)^(2) 

a A 

we obtain 

D[$]D[<7]D[g$]D[g$]D[g CT ]D[g CT ] 

exp 

i JV 

- « E (t-^! 2) + - 2) + 2)) $ a (2) 



(A.40) 



dld2 fyg$(l,2)g$(l,2) + i (iVg ff (l,2) - $> A (l)<7 A (2))g„(l,2) 



exp 



a = l 

1 



- |J1J2 ((2J o g* 2 (l,2) + «5(l-2)) g CT (l,2) 



(A.41) 



Making a shift (A. 34) over Q§ and the integration over $ and a: 
N 



exp 
exp 



2 

N 



j d\d2 (g $ (i,2)g $ (i,2) + g ff (i,2)g CT (i,2) + (- J Di 2) + /U (t)) g$(i,2)) 



TrLnQ $ + TrLnQ a 



N 



dld2 2J o g^(l,2) + £(l-2) g CT (l,2) 



2 v -v* -v V 2 

Eliminating Q $ and through saddle point evaluation, we obtain the action: 



(A.42) 



D[Q$]D[Q a ] exp(-NS) 
2S= f dld2 \{-D { ? + //(*))g$(l,2) + [2J g* 2 + 6] . g CT ] - TrLn[g $ ] - TrLn[g ff ] . (A.43) 



The equations of motion are obtained by saddle point over Q§ and Q a 

= (~d{ 2) + //(*)) 5(1 - 2) + 4J g$(l, 2)g CT (l, 2) - Qs\l, 2) 

g; 1 (i,2) = (2j g- 2 + «5)(i,2) 

and multiplying the first equation by Q§ and the second by Q a we finally get: 

o = (-d[ 2) + ^it)) g$(i, 2) + 4j ((g$ . g CT )g$)(i, 2) - s(i - 2) 
<5(i-2) = ((2j g* 2 + £)g ff )(i,2). 



(A.44) 
(A.45) 



(A.46) 
(A.47) 



These equations when written in components become equivalent to the equations of motion 
for the SCSA, Eqs.(III.4)-(III.8). 
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The generalization of what we have done to several modes that derive from an energy is 
straightforward. 

We can now see the formal difference between MCA and SCSA. Both the MCA and the 
SCSA equation for Q§ are of the form: 

{-D (2) + fi(t))Q(l, 2) = 8(1-2) + J d3 m[Q](l, 3)Q(3, 2) (A.48) 

The kernel m[Q] is a function m[Q](l } 3) = V"(Q(1,3)) for MCA, while it is a non-local 
functional m[Q$] = 4J (2J + <5) _1 • Q$ in SCSA. 

Symmetries 



Let us finally turn to the question of the symmetries associated with the equilibrium 
theorems. The SUSY group is generated by three operators 40-42 : 

D' = T§- e + e§- t 

= m (A-49) 



D' 2 = D' 2 = (A. 50) 

We can construct a version of this group that acts on two-point (in general n-point) functions, 
as: 



D' = D'(l ) + D'{2) 
D' = D'(l) + D'(2) 

[D',D'] + = ^ + ^ (A.51) 

The meaning of the three generators can be understood when they are made to act on a 
correlation function. Firstly, causality plus probability conservation imply (irrespective of 
equilibration): 

D'Q(f,2) = (A. 52) 

(This was already assumed in selecting the non-zero terms in (A. 10)). The other two gen- 
erators: 

(£ + £) w - 2 > = ° ~ TT '' 

D'Q(1,2) = -»■ FDT . 

The question of non-equilibrium can now be stated as follows: 

• Systems with explicit non-equilibrium have a dynamical action (or dynamical equations 
of motion) that break SUSY explicitly. 
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If the system is a priori able to achieve equilibrium, then SUSY is not explicitly broken. 
Then, two things may happen: 

• The effect of the initial conditions is a finite transient ~ t eq in which TTI and FDT do 
not hold. In this language, SUSY is unbroken by the boundary conditions. 

• If the system never achieves equilibrium, as in the case of the low-temperature version 
of the MCT equations t eq — > oo, the effect of the initial conditions is then to break SUSY 42 ' 45 
(violate TTI and FDT) well within the 'bulk' of times. SUSY is then spontaneously broken. 

The initial conditions play for SUSY (FDT and TTI) the same role played in ordinary 
symmetry-breaking by space boundary conditions: if the symmetry is spontaneously broken 
their effect extends away from them. 

Hence, if one is treating a system like a spin or structural glass within an approximation 
(or a phenomenological model), one must make sure that the resulting theory does not break 
SUSY explicitly, otherwise one may be introducing non-equilibrium by hand. 
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Appendix B: 



Derivation of the SCSA equations from a disordered Bernasconi model 



In this Appendix, we derive the dynamical equations corresponding to the disordered 
Bernasconi model using standard functional methods. For a more compact derivation using 
supersymmetric functional methods, see Appendix A. Following Ref. [11], we will define the 
disordered version of the Bernasconi model by the following Hamiltonian: 

T N N , 

u = EtE J«,f^M 2 + f E € + ^ E 

For each A independently, j£ @ = Jp a is equal to 1 with probability and zero otherwise. 
Each A thus corresponds to a certain 'pairing' of the 'sites' {a}. For J > 0, the ground 
state corresponds to a configuration of the <f> a which simultaneously minimizes all the 'partial 
correlation' S\ defined for each pairing as 

1 N 



S &) = -Iff E J i^h (B.2) 



N a<fj 

The problem is extremely frustrated if the <f> a are Ising spins, and becomes trivial if the <f> a 
are unconstrained continuous variables. The model defined in (B.l) is intermediate, since, 
as usual, a certain amount of constraint is enforced by the terms proportional to // and g' . 
The dynamical equations read: 

^ = ~ + Vc(t) (V a (t)vp(t')) = 2T6 a , p 6(t - t') (B.3) 

Now, the two identities (B.2),(B.3) can be written in a convenient functional way by 
expressing the 8 functions in Fourier space. (We set g' = for simplicity, and will discuss 
the modifications induced by g' ^ at the end of the calculation). We thus write: 

n[n^«(o#«(t)][n^A(t)^A(t)]exp- / AE^i^)— e ^v-wmoh 

t a A A V^v a <f3 

+ £ Ut){k(t) + nM*) + 4E s x (t)j^Mt) + ii*(t)}] = i (ba) 

V A Aj/3 

(The Jacobian associated with the change of variables is equal to 1 if one uses the Ito 
prescription). Averaging over the (Gaussian) thermal noise amounts to replace <f> a (t)r) a (t) 
by T(f) a (t) 2 . The average over the j£ ^ can also be performed and leads, for A" large, to the 
following expression: 

eX PEEt"^ / dt(-S\(j) a (j)i3 + JoSx^a^fj + JoS\<j) a <i){3){t) + 

A a</3 
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-^-^ J dtdt'(-S\(f) a (f)fj + J S\(f) a (f)fj + J S\(f) a (f)fj)(t)(-S\(f) a (f)fj + J S\(f) a (f)fj + JoS\(j) a ^i3)(t')} 

The second is a sum of iV 3 terms with a positive mean, and is thus of order iV, while the 
first term is a sum of iV 3 terms of random sign, and is of order 1, which we neglect. 

The next step is to define six 'correlation functions', associated to the fields S\, S\. 

Let us introduce: 

c^t, f) = n- 1 Mt)Mt') G^t, f) = n- 1 Mt)k(t') 

a a 

z ct (t,t') = N- 1 J2k(t)k(t') (B.Q) 

a 

and similarly for the Cs 7 G S} Zs- The expectation values of G are actually response 
functions 44 and those of Z are in fact zero 44 but Z must be kept in the intermediate steps 
of the calculation. Again, these identities are expressed as 8 functions, introducing six new 
'conjugate' variables NCs,Si NGs,s^ N Zs,s- The 'interaction term' (B.5), expressed in terms 
of the C, G, Z simply reads: 

exp / dtdt'[Z s (t, t')Ct(t, ff + 2J 2 C 5 (t, t')Z 6 (t } t')Ct(t, f)- 

2 Jt'<t 

2J G s (t , ,t)G^t,t , )C^t,t') - 2J Gs(t,t')G^t',t)C^t,t , )+ 

+2 JlG^t', t)Gt(t, t')C s (t, t')] (B.7) 

The point now is that all the terms cont are proportional to iV, and can be 

treated within a saddle point approximation which becomes exact when N is large. The 
saddle point equations read: 

d 

— > C<f, = ZsCtf, + JqZ^Cs — JqG\G& — JqGsG\ (B.8 — 



dC, 



d 



dGs 

d 



— > Gtf, — —JqG^Cs + JqG^Cs (B.8 — b 



dZ s 

d 



— > Zs = JnCsCj, (B.8 



dC s 

d 



— > Cs = JqG^Gs — JqZsCs (B.8 — d 



dG ^ --0^G s = -JoGlCs (B.8 -e 



d „ A c, 



2 



— > Z s = (B.8 



dZ s 2 

where we have dropped the arguments (t,t') when they appear in the correct order (t > t' 
and indicated with a f when they appear in reversed order. 
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From their physical interpretation (see, e.g. Ref. [44]), one expects that Zs = Z^ = 
and Gtf, = Gs = for t < t'. The saddle point estimate of Eq. (B.4) averaged over the 
J's then finally reads: 

/ ]\[]ld(f>4t)d^t)][]ldSx(t)dSx(t)]exp- [ dt[J2S x (t){S x (t) + 

J t a A J A 

/ dt'J C^Sx(t') - l -ClS x {t')} + k(t){k(t) + P4a{t)+ 
Jt'<t 2 

+ / dt'MC^Gs - JoCsG^M?) ~ (TS(t ~ t') + JlCtCs)^')}} = 1 (B.9) 
Jt'<t 

This last equation is easy to interpret by comparison with Eq. (B.4) - it is simply the 
Fourier representation of some 8 functions implementing the following equations of motion: 

Gzlm = J j*dt'(J CsG4, - C4GsW) + £(t) (5.10) 

with G~l = f t + fi and {£{t)£{t')) = 2T8(t - f) + J 2 C 5 C^ and 

GosS(t) = -Jo fdt'C^Sit') + ((t) (5.11) 
Jo 

with Gos = 8(t — f ) and (C(t)C(t')) = |C|. Note that C, G are self-consistently determined. 
Hence, comparing with Eqs. (III. 10, III. 11), we indeed see that the SCSA for the usual </> 4 
theory are the exact equations describing the random Bernasconi model defined by Eq. 
(B.l), provided one identifies J with the </> 4 coupling, The only missing part is the 
'tadpole' contribution, which can be easily added by choosing a suitable value of g 1 in Eq. 
(B.l), since this last term only adds in the equation of motion of <f> a a non fluctuating 
contribution —g'cf) a C^{t,t) (for N large). 

One should note that, as emphasized in the text, g > in the original </> 4 theory cor- 
responds to J > 0, and hence to a well defined (bounded from below) Hamiltonian 7i. 
The dynamical SCSA equations are thus expected to have sensible solutions for all values 
of parameters, contrarily to the direct MCA. 
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